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Disordered Fermi-Dirac distributions are used to model, within a straightforward and essen- 
tially phenomenological Boltzmann equation approach, the electron/hole transport across graphene 
puddles. We establish, with striking experimental support, a functional relationship between the 
graphene minimum conductivity, the mobility in the Boltzmann regime, and the steepness of the 
conductivity parabolic profile usually observed through gate-voltage scanning around the charge 
neutral point. 
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I. INTRODUCTION 

A great deal of effort has been currently devoted to 
the study of the transport properties of graphene, with 
particular focus on the so-called minimum conductivity 
problem, i.e., the fact that monolayer graphene's conduc- 
tivity drops to a finite minimum value when the Fermi 
level is set at the charge- neutral point, where the density 
of states should hypothetically vanish (for non-disordered 
and non- interacting electrons) (lH5|. A careful attempt 
to use standard linear response theory in this context is, 
however, plagued by ambiguities related to the alterna- 
tive prescriptions that can be used in the computation of 
quantum matrix elements @. Notwithstading such puz- 
zling theoretical issues, it is generally accepted, on clear 
experimental grounds, that the minimum conductivity, 
for both mono and bilayer graphene, is a non-universal 
quantity (frequently observed to be around 4e 2 //i), being 
essentially dependent on the amount and the nature of 
electronic disorder. 

It is a spread view that a solution of the minimum 
conductivity problem has to do more with finding ap- 
propriate models of disorder than with improving elec- 
tronic transport theory. In this paper, we rely on the 
"puddle picture" of disorder, which is known to yield 
a consistent explanation of the finite minimum conduc- 
tivity phenomenon for a large class of graphene systems 
[E ll "El I ■ The essential idea of the puddle picture is that 



the disorder produced by the sample and substract im- 
purities is associated to local and smooth shifts of the 
electron energy spectrum. Due to the uniqueness of the 
Fermi level, the electric charge is distributed in disor- 
dered graphene in the form of negatively (electrons) or 
positively (holes) charged puddles. Furthermore, it is be- 
lieved, in connection with the celebrated Klein tunneling 
effect [l3j , that electric charge can be transported across 
the puddles with strongly supressed backscaterring. At 
the charge neutral point, graphene is then depicted as a 
mixture of electrons and holes, a fact that was actually 
pointed out in the very first graphene transport experi- 
ments (l4| . 

We note, as an important remark, that the graphene 
conductivity has a non- vanishing minimum value even in 



clean suspended samples, where charged puddles should 
not exist due to the absence of a substract. Thus, the 
puddle picture cannot be evoked in these cases, where it 
is likely that more general quantum phe nomena play a 
relevant role in the transport process |15Hl7j ]. 

It is clear that a comprehensive test of the puddle pic- 
ture should address, with reasonable accuracy, the char- 
acterization of disorder in the graphene samples under 
investigation. Having in mind that this is a modeling 
task of difficult validation, it would be of great interest 
to establish general results that would not rely on most 
of the disorder details. This is our aim in this work, to 
be pursued within the framework of the Boltzmann ap- 
proach, however through an alternative implementation 
of the puddle picture of disorder. 

This paper is organized as follows. In Sec. II, we 
rephrase the puddle model of disordered graphene as 
a network of coupled fermion systems with randomly 
shifted energy spectra. This system is studied through 
the usual Boltzmann transport theory in Sec. Ill, where 
our main result - a statement on the behavior of con- 
ductivity near the charge neutral point - is derived and 
clearly confirmed from available experimental data. In 
Sec. IV we discuss, with the help of standard phenomeno- 
logical ideas, why the application of the Boltzmann for- 
malism is indeed meaningful in the analysis of the con- 
ductivity profile around the charge neutral point, a fact 
that could seem paradoxical at first sight, since near 
the Dirac point the Fermi wavelength is large enough 
to break the semiclassical transport regime. In Sec. V, 
we summarize our findings and point out directions of 
further research. 



II. STATISTICAL DESCRIPTION OF 
GRAPHENE PUDDLES 

Consider a system of gapless fermions with one-particle 
spectrum e^. = ak a (a — 1 and a = 2 correspond to ideal 
monolayer and bilayer graphene systems, respectively). 
As sketched in Fig. 1, we depict graphene electron- 
hole puddles as open ideal, two-fluid subsystems sharing 
chemical potential fj, for electrons and —\x for holes, in a 
large charge transport network 18]. The role of disor- 
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der is encoded in the "energy broadening function" 
which is just the probability density function of finding a 
subsystem with shifted electron or hole energy spectrum 
€k — £ in a given puddle. The energy broadening width 
Seo is defined from 



(1) 



and is assumed to yield a complete parametrization of 
In other words, we have 



de 



(2) 



where g(-) is a dimcnsionless universal probability density 
function (an educated guess is to take it as a zero-mean 
gaussian distribution with unit standard deviation, but 
its exact form is not important in our discussion). 




FIG. 1: Schematic view of disordered graphene as a network 
of open smaller ideal subsystems Qi (electron- hole puddles), 
which have one-particle energy spectra et — £», where the £j's 
are independent and identically distributed random variables. 
Arrows indicate the flow of particles and holes through the 
boundaries of f2i. 

In accordance with the statistical picture of electron- 
hole puddles advanced here, we put forward the zero- 
temperature equilibrium electron/hole occupation num- 
bers at wavenumber k as 



dep(oe(e±A*-e*) 



(3) 



where O(-) is the Heaviside function, and the above pos- 
itive and negative indices refer to electrons and holes, re- 
spectively. We have tacitly assumed, in Eq. ©, that 
fermions are described by an electron-hole symmetric 
hamiltonian [l9j . 

Observe that the total electric charge 



= e (^) / d 2 k[f+(k)-f-(k)] 



(4) 



vanishes, from Eq. ([3J, at /z = 0, the so-called charge 
neutral point. On the other hand, the carrier density 
(i.e., the mean number of electrons and holes per unit 
area) at the charge neutral point, no, is finite, and is 
readily evaluated as 

n = n + (0) + n-(0) = 



{2-Kf 
2tt \ a 



dfrKO J - ak a ) 



(5) 



III. CONDUCTIVITY PROFILE AROUND THE 
CHARGE NEUTRAL POINT 



Under the action of a small external electric field E, 
time-dependent occupation numbers can be obtained in 
principle as solutions of the usual Boltzmann equations 
in the relaxation time approximation, viz., 

[|±^-V fe ]/ ± (fc,t) = -^[/ ± (fc,t)-/ ± (fc)] , (6) 

where is the scattering time at wavenumber k. The 
stationary solutions of §6§ are, up to first order in the 
electric field, 



it 



eEr k 



(7) 



Substituting (0) in the expression of the electric current 
density, 

J= ih) 2 1 d2 k[f + ik)-r{kw-'k , (8) 



we find the conductivity 



h 2 



where 



+ m) + - M)]fo t , 



(9) 



(10) 



with fc s = (£/a) 1/Q . 

The scattering time at wavenumber k is expected to 
have the general form Tfe = F(k)/ni mp , where ni mp is the 
impurity concentration and F(k) is some scaling function 
of k 0, @|. Rather than attempting to compute F(k) on 
a first-principle basis, we determine how it should scale 
with k from the observed behavior of conductivity. More 
concretely, we assume that far from the charge neutral 
point, that is, for \/j,\/6eo ^> 1, the conductivity becomes 
a linear function of the carrier density n. In this region, 
(i) the chemical potential can be identified to the Fermi 
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energy of an ideal gas, that is, \fi\ — ak a , and (ii) Eqs. 
© and (TTU1) yield <r = onre 2 [iTk^/2h 2 . Since the carrier 
density is n oc k 2 , the conductivity will depend linearly 
on n only if F(k) oc k 2 ~ a . Therefore, we may conven- 
tionally write the scattering time as 



Tk 



chk 2 



(11) 



'imp 



where c is a dimensionless pref actor. It is worth of men- 
tioning that Eq. (TTT|) , which is taken to hold also around 
the charge neutral point, can be in fact derived for mono- 
layer graphene (a — 1) in the case of Coulomb impurity 
potentials [!, [|[. 

In order to investigate the conductivity behavior close 
to the charge neutral point, let us expand (flOl) around 
l-i = 0. We have, up to second order in ^i, 



t( m ) = 2^ oo ^[p(0 + yp"(0] 



(12) 



Substituting (fTT|) in (|12[) and taking ([2]) into account, we 
obtain 



2ch 



l imp 



so that 



where 



co 



9(0 



Sep 
a 

1 

h 2 



JL 

Se 



g"(0 



°± ( M \ 
2 \6e J 



CO 
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is the minimum conductivity, and 



o"i = era 



(13) 



(14) 



(15) 



(16) 



It is usual to control the two-dimensional carrier charge 
density in transport experiments through the voltage bias 
provided by gate electrode devices |l|. Variations of the 
gate voltage V g are proportional to charge density vari- 
ations CgAVg, where c g is the gate capacitance per unit 
area. Let us take, without loss of generality, V g = at 
the charge neutral point. We can relate /j, to V g from the 
expression of the total charge density e[n + (fj.) — rt_(/j)]. 
With the help of Eq. ([3]), we get, up to first order in ^, 



n+{p) - n-(fi) = 



1 



d£p(0 / d'k 



(2n) 2 

x[Q(Z + n-e k )-9(£-n-e k )] 

=ir*^- M) - p(e+/i)i (a 



+ o(m 2 ) 

de«?'(Oe- + ^(M 2 ) • (17) 
The above result implies, thus, that for <C Sep, 



4ira V a 



e/i / Sep 



(18) 
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FIG. 2: The dimensionless steepness, a" /a", of the graphene 
conductivity parabola near the charge neutral point is com- 
pared to the dimensionless combination of minimum conduc- 
tivity and mobility, aofJ.^ h /iJofi^, h . Open circles and crosses 
refer to electron and hole mobilities, respectively [20]. The 
other symbols are related to non-specified electron or hole 
mobilities. Data source: triangles [21| |; squares [%|; asterisks 
(bilayer graphene) [2^ ]. The squares and asterisks are slightly 
shifted to the right and to the left, respectively, to help visu- 
alization of the linear interpolation, which has unit slope as 
predicted by (|19[) . 



We have at this point all the ingredients to establish a 
phenomenological theorem for the steepness of the con- 
ductivity parabolic profile around V g — 0, viz., 



d 2 a 



oc aofJ-eh 



(19) 



v„=o 



where /i e /, = a /en oc l/nj mp is the carrier mobility mea- 
sured far from the charge neutral point. Relying now 
upon the fact that at the charge neutral point free charge 
carriers come from donor or acceptor impurities, we have 
7iQ OC Tlimp [1 01 ] . so that fx e h oc 1/no. Relation (|19j) . then, 
follows immediately from ©, ([Tg |) -([T5 ]) . and IfTg]), 

We have taken a", ao and /j, e h from the experimen- 
tal data reported in previous studies [20M23l j . For each 
group of measurements, carried out for various impurity 
concentrations, we denote by a" and ctq the parabola 
steepness and the minimum conductivity associated to 
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the largest selected mobility We work, then, with 

the dimensionless quantities a" /a" and o'ofJ-lh/^'oP'eh t° 
find, as it is shown in Fig. 2, compelling evidence for the 
validity of ([15]), 

It is interesting to check, additionally, if the statistical 
description of graphene transport that has led ultimately 
to (TT5)) is also able to provide a reasonable estimate of the 
minimum conductivity value. Note that with the help of 
Eqs. © and (fT5|), we find 

ao = 2nc^^-. (20) 

The dimensionless constant c is, in the case of Coulomb 
impurities, a function of the Wigner-Seitz radius r s , 



where 

G{r s )=r 2 ^~^[r 2 s g{2r s )]} (22) 

with 

7T 11 

g(x)=-l-{ t + (1 - t 2 ) , - arr.cos - (23) 

2 v^ 2 - 1 ^ 

Considering Si02 as the prototypical substract, we take, 
as estimates for monolayer graphene, r s — 0.8 H H[ 
and no = n imp /10 [l(|. It follows, from Eqs. ([20])- 
(|23| that these two numerical values conspire to give 
(To = 1.00834 x e 2 /h, which is incidentally very close to 
half of the conductance quantum. Taking into account 
the double valley and the spin 1/2 degrees of freedom of 
graphene electrons, its minimum conductivity is expected 
to be around 4e 2 /h, as in fact it has been reported in sev- 
eral transport experiments. To be fair, however, we stress 
that the measured CTo's are strongly sample-dependent 
[E HIIj varying typically in the range (2 — 7)e 2 /h. As 
we clarify in the next section, our evaluation of the mini- 
mum conductivity has to be interpreted more as an order 
of magnitude estimation than as a procedure to obtain a 
numerically precise result. 

IV. CRITICAL REMARKS 

We have found, in agreement with experiments, that 
the conductivity depends linearly on the density of charge 
carriers n for 3> Seo, where as for <C Seo, the 
conductivity profile has a parabolic shape. However, one 
may get puzzled by the fact that both (p~9|) and the above 
estimate of the minimum conductivity have been derived 
within the Boltzmann approach, while, as a matter of 
principle, the transport physics around the Dirac point 
is not semiclassical, but associated to relevant quantum 
corrections [24] . A second source of difficulty is that pre- 
vious semiclassical analysis, as the ones provided by the 



self-consistent Born approximation (25j and from the nu- 
merical diagonalization of disordered hamiltonians [26| . 
actually lead to minimum conductivities which are in 
general smaller than 4e 2 /h. 

In order to address a solution of this confusing state of 
affairs, we note, as a key point, that Seo can be identi- 
fied, within the puddle picture of disorder, with the en- 
ergy scale where the semiclassical-to-quantum crossover 
takes place. To understand it, first recall that the en- 
ergy broadening Seo is ultimately due to electrons deliv- 
ered or captured by the substract impurities; such an 
energy scale is thus related to the carrier density no, 
which is on its turn proportional to r^ mp , as already 
discussed at the end of Sec. III. The typical Fermi 
wavelength at the charge neutral point is, furthermore, 
Xf ~ 1/y/no ~ l-fy/nimp. In second place, we call at- 
tention to the fact that the mean distance between the 
scattering impurities is 8 ~ 1/ ^Jni mp ~ A/ , which indi- 
cates that Seo gives indeed a crossover energy scale for 
electronic transport in graphene. 

If we assume now that transport in the chemical po- 
tential range |<5/z| > Seo is to some good approximation 
described by the formalism developed in Sec. Ill, then it 
is clear that the analytical matching between the semi- 
classical (linear) and quantum (parabolic) profiles of the 
conductivity at the energy scale Seo is the essential reason 
for the validity of (flU| . From this point of view, we also 
conclude that our estimate of the minimum conductivity 
is not quite a direct semiclassical evaluation - rather, it is 
the result of an extrapolation, by curve matching, of the 
conductivity behavior at the crossover energy scale Seo- 



V. CONCLUSIONS 

We have introduced a statistical version of the pud- 
dle picture of graphene disorder, which provides a simple 
and straightforward phenomenological Boltzmann trans- 
port treatment of graphene conductivity. We have been 
able to predict and verify in this way an unsuspected re- 
lationship that holds for the geometrical parameters that 
define the parabolic shape of the conductivity profile near 
its minimum value, as a function of the backgate voltage, 
and the mobility of charge carriers in the linear regime, 
far from the charge neutral point. 

The experimental validation of the original formula- 
tion of the puddle model 0] is an involved issue, which 
goes beyond the mere fitting of measured conductivity 
profiles. Our central result, Eq. ([!§]) , bypasses this kind 
of difficulty, since it does not rely, in principle, on specific 
modeling details, like the density, statistical distribution 
and nature of the disordering impurities. Fig. 2 provides, 
therefore, an important support for the general validity 
of the puddle picture of graphene disorder. 

Natural extensions of the present work are related to 
the consideration of electron-hole asymmetry, tempera- 
ture, and interaction effects. It would be interesting to 
investigate if similar results can be derived for the case 
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of clean suspended graphene samples, where the puddle 
picture of disorder is probably not adequate anymore. 
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